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We describe a procedure for mapping a self-consistent mean-field theory (also known as den- 
sity functional theory) into a shell model Hamiltonian that includes quadrupole-quadrupole and 
monopole pairing interactions in a truncated space. We test our method in the deformed N — Z 
sd-shell nuclei 20 Ne, 24 Mg and 36 Ar, starting from the Hartree-Fock plus BCS approximation of the 
USD shell model interaction. A similar procedure is then followed using the SLy4 Skyrme energy 
density functional in the particle-hole channel plus a zero-range density-dependent force in the pair- 
O j' ing channel. Using the ground-state solution of this density functional theory at the Hartree-Fock 

D ' plus BCS level, an effective shell model Hamiltonian is constructed. We use this mapped Hamilto- 

00 , nian to extract quadrupolar and pairing correlation energies beyond the mean field approximation. 

The rescaling of the mass quadrupole operator in the truncated shell model space is found to be 
almost independent of the coupling strength used in the pairing channel of the underlying mean-field 
theory. 
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PACS numbers: 21.60.Jz, 21.60.Cs, 21.10.Dr, 21.30.Fe 



I. INTRODUCTION 



One of the major challenges in nuclear many-body theory is to develop a microscopic and systematic approach that 
would account for the most relevant correlations in the description of nuclear properties. Most existing methods can 
be divided into two major classes: the self-consistent mean field (SCMF) approximation also known as density 
functional theory (DFT), and the configuration- interaction shell model (C1SM) approach Q. SCMF theories are often 
■ cast in terms of an energy density functional that is minimized to obtain the ground-state solution of the system. 
There exist parameterizations of this energy density functional that are valid globally through the table of nuclei. 
Q\ • Such parametrizations are usually based on the zero-range Skyrme force Q or on the Gogny interaction Q • The CISM 
approach requires as input a large number of interaction matrix elements in a truncated space. While effective CISM 
interactions can sometimes be traced back to the bare nucleon-nucleon interaction, it is often necessary to adjust them 
empirically. Since the CISM method requires a truncation to a finite number of shells, the single-particle energies and 
^ , interaction matrix elements are specific to the mass region under consideration. Thus the CISM approach lacks the 
global validity of the SCMF. 

On the other hand, the CISM approach has the advantage that the underlying nuclear wave functions (both ground 
and excited states) are fully correlated. Complete Ofrcu major shell calculations in p- d, sd- dU and /p-shell @,H 
nuclei have successfully described many of the observed properties of these nuclei. In contrast, the SCMF approach 
does not account for all correlations present in the nuclear ground state. Partial correlations in a mean-field approach 
are often accounted for by the spontaneous symmetry breaking of the ground state (e.g., the breaking of rotational 
symmetry by a deformed ground state) 0] . Additional correlations beyond the mean- field approximation are included 
by restoration of this broken symmetry (e.g., by angular momentum projection) l lOf , and by configuration mixing 
of symmetry-projected states via the generator coordinate method (GCM) [ll], Il2l. Il3j. The GCM is also useful in 
providing information on excited states of the nucleus. Angular-momentum and particle-number projected GCM has 
been applied successfully in global studies of binding energies [ljj [l6[ and of spectroscopic properties of the first 
excited 2 + state [ItJ in even-even nuclei. 

In a recent work pj| we initiated an approach that takes advantage of the global validity of the SCMF method 
and the higher accuracy of the CISM approach. It is based on the idea of mapping an SCMF theory onto an effective 
CISM Hamiltonian that is defined in a suitably chosen truncated model space. Such an approach would enable us 
to start from an SCMF theory with a global parametrization and construct effective CISM Hamiltonians in different 
mass regions. In heavy nuclei, the required truncated shell model spaces are expected to be too large for conventional 
diagonalization of the CISM Hamiltonian. In such cases, one can apply methods that have been developed to treat very 
large model spaces, such as the shell model Monte Carlo (SMMC) method [l|J[2(|j * ne Monte Carlo shell model [21| 
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approach, and direct matrix methods [22j |. Alternatively, one can use tractable approximations that go beyond the 
mean- field approximation such as the RPA [23| . 

A mean-field approximation, based on a Skyrme force, was used in Ref. to determine the single-particle energies 
of a CISM Hamiltonian. In Ref. [ID we carried out the first step of mapping an SCMF theory onto the CISM and 
constructed an effective Hamiltonian that includes quadrupolar correlations. The map was tested successfully in 
N = Z s<i-shell nuclei, using the Hartree-Fock (HF) approximation of the so-called USD Hamiltonian [6| as our 
starting point for an SCMF theory. The USD Hamiltonian was constructed empirically and is known to give a 
good description of the low-energy spectroscopy of se?-shell nuclei. Using the mapped CISM Hamiltonian, we found 
that quadrupolar correlations alone account for more than 50% of the full correlation energies contained in the USD 
Hamiltonian. We then applied our mapping procedure to study quadrupolar correlations in N = Z sd-shell nuclei 
starting from the SLy4 parametrization of the Skyrme energy functional as the SCMF theory [l8[ . 

Here we explore the SCMF to CISM mapping idea further by including pairing correlations in the mapped effective 
CISM Hamiltonian in addition to quadrupolar correlations. The starting point is an SCMF theory that includes pairing 
correlations. Such a theory can be solved in the HF plus BCS approximation [or more generally in the Hartree-Fock- 
Bogolyubov (HFB) approach] to provide the ground-state single-particle HF wave functions and energies, as well as 
the BCS occupation amplitudes. A similar HF+BCS approximation can be used to determine a spherical single- 
particle basis. This mean-field output is used to construct an effective CISM Hamiltonian in a truncated spherical 
model space that contains both quadrupolar and monopole pairing correlations. After testing the method using 
the HF+BCS approximation of the USD interaction as our initial SCMF theory, we apply it to deformed N = Z 
sd-shell nuclei in the context of DFT. The latter is based on the SLy4 parameterization of the Skyrme force in the 
particle- hole channel [25| plus a zero-range density-dependent force in thepairing channel (26[. It is solved in the 
HF+BCS approximation using the Brussels-Paris computer code ev8 [13, [H| . The mapping constructed here could 
also be useful in global studies of level densities within the SMMC approach [29|, [3(| ■ Both quadrupolar and pairing 
correlations play an important role in the microscopic calculations of level densities. 

Pairing correlations are natural for spherical nuclei. Here, however, we consider the addition of pairing correlations 
on top of quadrupolar correlations in the CISM Hamiltonian. We therefore focus on nuclei whose HF+BCS ground 
state is deformed. In principle, the method presented here can be extended to spherical nuclei by applying an external 
quadrupole field, thus effectively deforming the nucleus. 

The outline of this paper is as follows. In Sec. |TT] we discuss our mapping procedure when pairing correlations are 
included in the SCMF theory. The ground state is solved in the HF+BCS approximation and its unique description 
requires both a density matrix and a pairing tensor (or an anomalous density matrix). The CISM Hamiltonian is 
constructed in a truncated spherical space to include both quadrupolar and pairing correlations. In Sec. Mil we test 
our mapping procedure using the HF+BCS approximation of the USD Hamiltonian as our initial SCMF theory for 
deformed N = Z sd-shell nuclei ( 20 Ne, 24 Mg and 36 Ar). In Sec.[TV]we map a DFT that is based on the SLy4 Skyrme 
interaction plus a zero-range density-dependent pairing force onto an effective CISM Hamiltonian for deformed N = Z 
sd shell nuclei. The correlation energies extracted from the mapped Hamiltonian compare well with the correlation 
energies obtained with other methods. We conclude with a brief discussion of possible future directions in Sec. [V] 



II. THEORETICAL FRAMEWORK 



In this section, we discuss the theoretical framework of mapping an SCMF theory into a CISM theory in the 
case where the SCMF includes pairing correlations. We assume that the SCMF theory is solved in the HF+BCS 
approximation 9] . This is for example the case of Skyrme-type density functional theories in which pairing correlations 
are considered through the addition of a zero-range density-dependent force in the pairing channel (27l. |28|. 

The HF+BCS mean-field theory provides us with a set of single-particle HF orbitals 4>k and their corresponding 
single-particle energies and BCS occupation amplitudes Vf., as well as the total energy 

E mi +BCS of the HF+BCS 

ground state. For simplicity we assume that this solution, obtained as the global minimum of an energy density 
functional, is deformed. 1 Using the same SCMF theory (i.e., HF+BCS) we can also find the spherical solution. Our 
plan is to use this output of the SCMF theory to construct an effective CISM Hamiltonian that is defined within 
a truncated shell model space. In particular, we will consider an CISM Hamiltonian that contains a quadrupole- 
quadrupole interaction and a monopole pairing interaction. 

The CISM is defined in a spherical basis, and we determine such a basis from the spherical solution of the SCMF 
theory. A spherical solution could be found by iterations when the initial guess for the density matrix is chosen to 



1 Spherical nuclei can be treated by adding a constraining quadrupole field, leading effectively to a deformed ground state. 
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be spherically symmetric, and as long as the spherical symmetry is preserved in the iteration procedure. In Ref. [l8| 
we used the HF approximation, in which case it was necessary (for an open shell nucleus) to use the uniform filling 
approximation or add a pairing-like interaction with a fixed gap to preserve spherical symmetry in the next iteration. 
Here, the BCS approximation guarantees that the density matrix in the next iteration remains spherical. 

The construction of a standard spherical basis with good angular momentum j and projection m was discussed in 
Sec. IV of Ref. Il8l Here we use a as a generic label for the spherical basis, i.e., a includes the quantum numbers j, m 
as well as other quantum numbers. 2 

The transformation matrix between the deformed single-particle basis k and the spherical single-particle basis a is 
denoted by Uka = (k\a). We have 

where a J, is the creation operator in a spherical state a and b\. is the creation operator in a deformed state k. 



A. Deformed SCMF solution in the HF+BCS theory 

The deformed ground state in the HF approximation is a Slater determinant n^al]— ) which is uniquely described 
by a density matrix p. However, in the HF+BCS approximation, the ground state $ = Hk>o(uk + v ka\a^)\— ) is 
defined by the BCS amplitudes Uk,Vk (the product runs over half the number of single-particle states k > and 
k is the state conjugate to k). The unique description of this BCS ground state requires both a density matrix p 
and a pairing tensor (or anomalous density) k In the deformed basis ptk' = (SLSfc) and Kkk' = (bk'bk), where 
(. . .) denotes the expectation value in the HF+BCS ground state In the spherical basis, we can similarly define 
Paa' = (o-a'^a) and K aa / — (a a ia a ). 3 Using the unitary transformation (p} from the deformed to the spherical basis, 
we have 



= Yl Uh ' a ' pkk ' u kot > ( 2a ) 

kk' 

= ^ta > ( 2b ) 



kk' 



or, equivalently, in a matrix notation 



p -> U f pU , (3a) 
U^U* . (3b) 



K 



According to |3]), p transforms as a linear operator under the unitary transformation U while k transforms as an 
antilinear operator. 

Just as in the HF approximation, the trace of p must be equal to the total number of particles A of a given type 
(protons or neutrons) 

tr p = N . (4) 

In the framework of the HF+BCS approximation, this condition is satisfied by using a suitable chemical potential /i 
(for each type of particle). 

In the HF approximation, the ground state is a Slater determinant and the corresponding density matrix must be 
a projector p 2 — p with all eigenvalues cither 1 or 0. However, in the HF+BCS ground state the matrix p is no 
longer a projector and its eigenvalues v 2 , are between and 1. The matrices p and k, which characterize uniquely the 
HF+BCS ground state, satisfy now the relations 

p p 2 = kk) ; (5a) 
pK = Kp* . (5b) 



2 In general we use lowercase greek letters a, a' , ... to denote spherical states and lowercase roman letters k, k' , . . . to denote deformed 
states. 

3 For simplicity of notation, we use the same symbols to denote the matrices in the deformed and spherical basis and they can be 
distinguished according to their subscripts. 
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In the following, we assume that the SCMF theory has time-reversal symmetry so that the single-particle orbitals 
come in degenerate time-reversed pairs |fc) and \k) = T\k) (T is the time-reversal operator). In the following we choose 
the k > orbitals to be the orbitals with positive z-signature [2a ]. The set {k,k} with k > spans the complete 
single-particle space. We also assume that the ground state solution is axially symmetric, in which case J z — m is a 
good quantum number (z is the symmetry axis), and the positive z-signature orbitals have m = 1/2, —3/2, 5/2, . . .. 
Similarly, the spherical states appear in degenerate time-reversed pairs {a, a = Ta} with a > denoting the positive 
z-signature states. 

The matrices p and k have the following simple form in the deformed HF basis (for k, k! > 0) 

Pkk> = Phk' = v k$kk> , Pkk> = Pkk> = ; (6a) 
K kk< = ~ K kk' = u k v k S k k' , K kk , = K Vk , = . (6b) 

The transformation matrix U is real (and thus orthogonal) and the matrices p and k remain real in the spherical 
basis. Since the spherical states also have good z-signature, the transformation U does not mix states with different 
z-signature and for a > we have \a) — X)fe>o U ka \k). Applying the time-reversal operator to this relation and using 
the fact that U is real, we find that the time-reversed states transform in exactly the same way |a) = J2 k >o U ka \k). 
We therefore have 

U ka = U k& = ; U ka = U M . (7) 
It follows that the density and anomalous density matrices have the block form 

P=( o A ) 5 K= (°B o) ' ^ 

where A and B are real matrices whose dimension is half the total number of orbitals. It follows from Eqs. and 
© that (for a, a' > 0) 

A aa > = ^U ka >vlU ka , (9a) 

fc>0 

B aa ' = ^2 U ka 'U k V k U ka ■ (9b) 
fc>0 

The matrix p is real symmetric while the matrix k is real antisymmetric. The form ([8]) is equivalent to 

Paa' — Paa' — ^aa' t Paa' Paa' ? (10a) 
^aa' — ^aa' — ^qq' : ^aa' — ^aa' • (10b) 

It is convenient to define the density and anomalous density operators by 

p = EO fc >^ fc i + i^w) ; ( lla ) 



fc>0 



\k)u k v k (k\ - \k)u k v k (k\) , (lib) 



such that their matrix representation in the corresponding basis gives us the matrices p and k 1 i.e., p aa ' — {(ApW) 
and Kaa> = (a\k\a'). Equations ([9]) and (fT0|) follow then immediately from the operator form (|11| . 

The basic relations ((5j) defining an HF+BCS state can be rewritten in terms of the matrices A and B as 

A - A 2 = B 2 ; (12a) 
AB = BA . (12b) 

A more compact way of representing the relatio ns fl5t or (I12p is to introduce a generalized density matrix 7L whose 
dimension is twice the total number of orbitals 0, [31l| 



71 = 



( a a' a a) ( a a>a a ) 



P 

-K* 



K 

1-p* 



A 








B 





A 


B 








-B 


I A 





B 








I A 



(13) 
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where I is the identity matrix in a space whose dimension is half the number of orbitals. Relations ([5]) or (|12p are 
equivalent to the condition that is a projector 

K 2 = 71 . (14) 

Thus the eigenvalues of Tt must be 1 or 0. This condition is formally similar to the one satisfied by p in the HF 
approximation except that there are now twice as many eigenvalues. In our case, the matrix 7L is real symmetric and 
its eigenvalues are pairwise degenerate. 



B. Criteria for choosing the truncated CISM space and the rescaling of one-body observables 

As long as the spherical basis a is complete, U is unitary and we have just generated a different representation of 
the same ground-state density matrix and pairing tensor. However, in constructing a shell model space, it is necessary 
to truncate to a valence subspace. We define such a subspace in terms of a projector P = ^2' a |a)(a| where Yl denotes 
a sum over a subset of the single-particle orbitals a. The density operator p and anomalous density operator k are 
then replaced by projected densities 

p -> PpP ; k -> PkP . (15) 

The corresponding matrices in the truncated spherical space are PpP and PkP. 

In choosing a suitable truncated subspace, we require the projected densities to preserve approximately their basic 
properties in the complete space. In the HF+BCS approximation, the trace of the density matrix is constraint to 
give the total number of particles [see Eq. and we require that a similar relation is satisfied approximately by 
the projected density PpP with the number of particles replaced by the valence number of particles 

tr(PpP) « Valence • (16) 

Another requirement is that the conditions ([5]) are satisfied approximately by the projected densities PpP and PkP. 
Alternatively, we require the matrices PAP and PBP to approximately satisfy (TT2]) . It is convenient to define the 
projected Tt matrix 

pnp = ( ppp PkP ^ (17) 

rwr ^ PkP P PpP J ' {U> 

and require that it satisfies approximately the condition (|14l) . namely that its eigenvalues are close to either 1 or 0. 
Formally, this condition is similar to the one we require from the projected density matrix p at the HF level [l8j |. 
except that the projected TL matrix is defined in twice the number of dimensions. 

In general, we can calculate the expectation value of a one-body observable O = ^2 aa i O aa ia^ a a a i in the truncated 
CISM space using the projected density 

{&) = Tt{6PpP) = ]T O aa , Pa , a = tr (OPpP) , (18) 

where the sum is restricted to the truncated spherical space. The ratio between the expectation value of O in the full 
and truncated spaces provides us with simple rescaling factors of one-body operators in the truncated CISM space. 



C. Effective CISM Hamiltonian 



In this section we construct the CISM Hamiltonian given the HF+BCS ground-state solution. In the self-consistent 
basis k of the HF+BCS, the HF Hamiltonian isAnf kk' ~ e ^kk' > where £fc are the energies of the deformed HF single- 
particle orbitals. We then proceed as in Ref. 18. The HF Hamiltonian matrix h m f is transformed to the spherical 
space and projected on the truncated CISM space. The projected one-body Hamiltonian is 

h = PU T h mf UP . (19) 

Next, we expand h in multipoles h = ~^2 K and construct the second quantized one-body tensor operators 

of rank K [see Eqs. (10) and (8) in Ref. [l8[. The lowest multipole is the monople — J2 a e « a a a ai which we use 
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to define the spherical single-particle Hamiltonian, while the K = 2 multipole is a quadrupole operator h^ 2 \ which 
we use to define a quadrupole-quadrupole interaction. 4 

We now consider the following effective CISM Hamiltonian 

H = m \g Q : ■ : -g P (P}P p + P,£P n ) (20) 

where :: denotes normal ordering and P^ is the monopole pair creation operator 

pt = £a+at, (2i) 

defined separately for protons and neutrons. Here gQ and gp are, respectively, quadrupole and pairing coupling 
constants. In comparison with the CISM effective Hamiltonian considered in Ref. [l8|, Eq. (|20| contains an additional 
monopole pairing interaction. We note that one might consider an isovector monopole pairing (which contains also a 
proton-neutron component). However, in the DFT applications we discuss in Sec. IIV1 the resulting HF Hamiltonian 
does not conserve isospin symmetry (because of the Coulomb interaction) and we only consider here an Hamiltonian 
of the form (|20[) . 5 Both quadrupole-quadrupole and pairing interactions in (|20[) are attractive (i.e., gQ,gp > 0), 
leading to an CISM Hamiltonian that has a good Monte Carlo sign in the SMMC approach [29| . 

The method outlined above determines the effective Hamiltonian (|20|) up to the two coupling constants gQ and gp. 
In Ref. [TH we determine the coupling constant of the quadrupolar interaction by matching the deformation energy 
of the SCMF theory with the deformation energy of the CISM effective Hamiltonian when the latter is solved in the 
HF approximation. When the effective Hamiltonian has the form (f20| . it is necessary to match two energy scales in 
order to determine the two coupling constants gQ and gp. One energy scale is the deformation energy, here redefined 
in the HF+BCS approximation 

^ +BCS = ^ +BCS ~ ^ +BCS , (22) 

where S™ +BCS is the energy of the spherical HF+BCS solution and E^J +BCS is the HF+BCS energy of the deformed 
ground-state solution. The deformation energy in (|22[) is usually smaller than the deformation energy calculated in 
the HF approximation. 

To introduce a second energy scale, we define -E pa ir to be the average interaction in the pairing channel 



£ pair = -rtr(«A) , (23) 



where A is the pairing gap matrix. 6 The energy -E pa i r can be calculated in both the spherical and deformed solutions 
and we define the pairing energy excess to be 

cpHF+BCS _ pHF+BCS _ pHF+BCS / 9zL \ 
0£/ pair — -^pair.mf "pair.sph ■ \ L ^> 

Since E pa i T is lower (i.e., more negative) in the spherical solution, the energy difference in (|24l) is usually positive. 

The two coupling parameters gQ and gp are determined by matching simultaneously the two energy scales (|22p and 
([24]) between the SCMF theory and the CISM Hamiltonian when the CISM is treated in the HF+BCS approximation. 

In general, the coupling parameters should be renormalized depending on the truncation used for the CISM space. 
These renormalization effects are included implicitly in our procedure by matching the energy scales calculated in the 
SCMF approach in the complete space with similar energy scales calculated in the truncated CISM space within the 
HF+BCS approximation. 

An important application of our mapping is in the calculation of correlation energies. The correlation energy E COTT 
is defined as the difference between the ground-state energy in the HF+BCS approximation to the CISM Hamiltonian 
and the fully correlated ground-state energy E gs of the CISM 

-Ecorr = -^mf +BCS ~ -^gs ■ (25) 



4 Our studies here are limited to even-even nuclei, for which time-reversal invariance implies even values of K. 

5 For the USD tests in Sec. IIIII we found that an isospin-invariant interaction constructed from an isovector monopole pair operator gives 
only a small correction to the correlation energy when compared with a Hamiltonian of the form H20I I. 

6 The gap matrix is defined by A;;/ = = ^Zfcfc' v u> ky K kk'> where v A is the anti-symmetrized two-body interaction. 
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Quadrupolar and pairing correlations compete with each other; the quadrupole-quadrupole interaction favors de- 
formation in which pairing is suppressed, while pairing favors the spherical configuration. If we consider the spherical 
configuration alone, then the BCS approximation already takes into account some of the gain in pairing correlation 
energy. However, in the deformed nuclei considered here we effectively gain more in mean-field energy by deformation. 
In the deformed ground state, BCS correlations are often significantly smaller than in the spherical solution (or even 
vanish). Nevertheless, The correlation energy described by (125|) is expected to increase when a pairing term is added 
to the quadrupolar CISM Hamiltonian. This increase arises from fluctuations of the pairing fields beyond the BCS 
approximation. In particular, these fluctuations sample the spherical configuration in which pairing correlations are 
significant. 

Even in spherical nuclei (which we do not treat here) we expect to gain a finite correlation energy in (|25p when 
adding a pairing term to the CISM Hamiltonian, i.e., there are pairing correlations beyond those described by the 
spherical BCS solution. The BCS approximation is valid in the limit when the pairing gap is much larger than the 
mean single-particle level spacing, and studies of pairing effects in ultra-small metallic grains revealed deviations from 
BCS theory in the crossover to the fluctuation-dominated regi me ]32l l33j . Nuclei belong to this crossover regime, in 
which the pairing gap is comparable to the mean-level spacing [3J]. 

In the remaining part of this section we further motivate our choice of the second energy scale (1241) . As already 
noted, pairing effects are weaker in a deformed configuration and could vanish altogether (i.e., k = 0) in the deformed 
ground state. For example, in our HF+BCS studies of the USD interaction in Sec. MI A[ we observe such a collapse 
of the BCS solution in all three deformed N — Z sd-shell nuclei. We then have 6E^+ BCS = -E*?+™? s , and only 

pan pan ,spn ' " 

the spherical solution is relevant for determining the second energy scale. An alternative second energy scale that 
could be used in such cases is — -£^ph +BCS ' the decrease in spherical energy in the HF+BCS approximation as 
compared with the HF approximation. In general 



pHF t^HF+BCS pHF+BCS 



S(h^)-lstv(p sph v A p sph ). (26) 



^sph sph pair,sph \ / 2 

where p sph is the spherical one-body density and v A is the anti-symmetrized two-body interaction. The gain in pairing 
energy from — E^^^h > is thus canceled in part by the increase S(h^) of the one-body energy and the increase 
^Str(p sph v A p spil ) of the average interaction in the particle-hole channel. Despite this partial cancelation, — E p ^^ s 
remains well correlated with E^ 1 — -£^ph +BCS an< ^' as l° n § as the deformed solution is unpaired, we find that both 
energy scales (|24p and (|26p lead to CISM Hamiltonians with similar coupling constants. 

However, the deformed ground state can generally have a non-vanishing (though weaker) pairing solution, as we 
observe in our DFT studies of Sec. IIV Al In such cases, we find that matching the energy scale (|26|) results in a 
correlation energy (f2l)]) of the CISM ground state that is unphysically large. We therefore adopt as the second 
energy scale to be matched when the deformed solution is paired. The scale (|2"4")l also has the obvious advantage that 
it is defined within a single SCMF theory (i.e., HF+BCS) and does not require the use of two different mean-field 
theories (i.e., HF and HF+BCS). 



III. TESTS FOR THE USD SHELL MODEL HAMILTONIAN 



In this Section we test the mapping procedure discussed in Sec.|H]by considering the HF+BCS solution of the USD 
Hamiltonian Q as our starting point of an SCMF theory. The USD interaction provides a rather good description of 
binding energies, spectra and electromagnetic transition intensities of sd-shell nuclei 0. In these studies we compare 
the correlation energies found from the mapped CISM Hamiltonian with the full correlation energies of the USD 
Hamiltonian. In Sec. IIV1 we apply the mapping formalism to a DFT that is defined by an SLy4 Skyrme force [25[ in 
the particle- hole channel plus a zero-range density-dependent force in the pairing channel [26j . 

Here and in Sec. [IV] we focus on N = Z nuclei, in which isovector quadruple-quadrupole correlations could be 
ignored. We also restrict our studies to nuclei whose self-consistent HF+BCS ground state is deformed. In the sd 
shell, this excludes 32 S whose self-consistent HF+BCS ground state is spherical. Spherical nuclei can be studied by a 
similar method once a deformation is induced by a quadrupole constraining field. 

We also restrict our studies to nuclei in which the spherical configuration supports a BCS solution (so that 
^paIr +BCS 7^ 0)- Otherwise, it is necessary to add a constraining pairing field. In the sd shell, this excludes 28 Si for 
which both spherical and deformed solutions are found to have no pairing solution. A similar situation occurs in the 
framework of the DFT application discussed in Sec. IIVI Therefore, the N — Z sd-shell nuclei we consider here and 
in Sec. |IV] are 20 Ne, 24 Mg and 36 Ar. 
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A. Mapping the HF+BCS theory of the USD interaction 



We have carried out HF+BCS calculations using the USD interaction. The resulting HF+BCS energies are listed in 
Table[T]for both spherical and deformed ground-state configurations. We find the ground-state solutions to be prolate 
for 20 Ne and 24 Mg, and oblate for 36 Ar, with average mass quadrupole moments (Q) given in Table HI The HF+BCS 
deformation energies (|22|) are calculated to be 5.64,7.05 and 1.27 MeV for 20 Ne, 24 Mg and 36 Ar, respectively. It is 
interesting to compare these energies with the corresponding HF deformation energies [l8[ of 10.92, 11.97 and 3.82 
MeV for the same nuclei. The reduction in deformation energies originates in a negative -E^ r + S ph S that lowers the 
spherical HF energy. 



TABLE I: HF+BCS results for the deformed N = Z sd-shell nuclei 20 Ne, 24 Mg and 36 Ar, using the USD interaction and the 
effective CISM Hamiltonian (|20p . Energies are in MeV, the average quadrupole moment (Q) is in b 2 (with b being the oscillator 
radius), gq is in MeV/b 4 and gp is in MeV. 



Nucleus 


Interaction 


9Q 


gp 


pHF+BCS 
^sph 


pHF+BCS 
pair.sph 


pHF+BCS 


pHF+BCS 
pair.mf 


rnHF+BCS 
-^def 


rpHF+BCS 
pair 


(Q) 


-C/corr 


20 Ne 


USD 






-30.74 


-7.04 


-36.38 





5.64 


7.04 


15.4 


4.1 




pairing 


0.0498 


0.2865 


-36.96 


-7.04 


-42.60 





fit 


fit 


14.9 


4.1 


24 Mg 


USD 






-73.12 


-7.00 


-80.17 





7.05 


7.00 


18.0 


6.9 




hr- 2 ' + pairing 


0.0254 


0.4765 


-103.42 


-7.00 


-110.47 





fit 


fit 


17.5 


5.7 


36 Ar 


USD 






-225.29 


-3.55 


-226.56 





1.27 


3.55 


-13.5 


4.0 




+ pairing 


0.0582 


0.3554 


-379.57 


-3.55 


-380.84 





fit 


fit 


-12.2 


2.7 



This leads us to the second energy scale SEflf BCS [Eq. f2"g|)]. which is found to be 7.04,7.00 and 3.55 MeV for 
20 Ne, 24 Mg and 36 Ar, respectively. The deformed solution of all three nuclei is found to have no pairing and, as noted 
in Sec. Ill C[ an alternative second energy scale is given by ([26]). 

The original space of the USD interaction is already truncated, so orbital truncation effects are not present. The 
transformation matrix from the deformed to the spherical basis is therefore unitary and the projector P is not 
required. In particular, the consistency relations (|16p and (O [or equivalently (|14[) ] are automatically satisfied within 
the spherical sd shell. 

The deformed mean-field Hamiltonian is expressed in the spherical basis [Eq. (|19p ] and decomposed into multipoles 
to find /i (0) and /i (2) . To construct the effective CISM Hamiltonian ([2D]), it remains to determine the strength 
parameters <?q and gp of the quadrupole and pairing interactions. We tune these parameters to match simultaneously 
both energy scales (|22|) and (|24| of the USD Hamiltonian when solved in the HF+BCS approximation. For example, 
in 20 Ne the USD energy scales of £ , JJ +BCS = 5.64 MeV and 5E^ BGS =7.04 MeV are matched by similar energies 
calculated for a CISM Hamiltonian (gD]) with g Q = 0.0498 MeV/6 4 and g P = 0.2865 MeV. Since the deformed 
configuration has no pairing solution, we can alternatively match the energy scale (|26j) (together with the deformation 
energy), resulting in very similar values for the coupling constants gQ and gp. 

To verify the validity of our map, we compare the quadrupole moments of the mapped CISM Hamiltonian with 
their mean-field values in the USD interaction. The results, shown in Table HI are in close agreement with each other. 

The ground-state energies of the CISM Hamiltonian (|2"0|) were calculated using the computer code oxbash [35[ and 
the corresponding correlation energies are listed in Table |U Figure Q] shows these correlation energies (solid squares) 
versus mass number A. They are compared with the correlation energies extracted from a mapped CISM Hamiltonian 
that includes only quadrupolar correlations (open squares) [TJ| and with the full USD correlation energies (open 
circles) . Quadrupole correlations alone seem to be responsible for more than 50% of the USD correlation energies [l8| • 
We observe that the inclusion of pairing correlations in the CISM Hamiltonian (on top of quadrupole correlations) 
brings us even closer to the full USD correlation energies. These results are encouraging considering that the USD 
interaction [6j with 3 single-particle energies 63 anti-symmmetrized interaction matrix elements was fitted to reproduce 
the observed properties of nuclei in this mass region. 
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^4 

o 

UJ°2 


20 24 28 32 36 
A 

FIG. 1: Correlation energies -E corr versus mass number A for the deformed N — Z sd-shell nuclei 20 Ne, 24 Mg and 36 Ar. Results 
for the mapped CISM Hamiltonian (120[) (solid squares) are compared with the quadrupolar correlation energies calculated in 
Ref. [HI] (open squares) and with the full correlation energies of the USD Hamiltonian (open circles). 



IV. MAPPING THE DFT 

In this section we apply the mapping procedure of Sec. [ill starting from a DFT that is based on the SLy4 Skyrme 
force [25[ in the particle-hole channel plus a zero-range density-dependent force 

V(r lt r 2 ) = -, 9 DFT (l - (l - ^) - r 2 ) (27) 

in the pairing channel In Eq. (f27|) . P a is the spin exchange operator and p(r) the total nuclear density. Through- 
out this work, we choose the strength of the force <? DFT to be the same for protons and neutrons and the central 
density is ,o c =0.16 fm -3 . A smooth cut-off of 5 MeV around the Fermi level is used for the interaction l|27p. The 
above DFT is solved in the HF+BCS approximation using the Brussels-Paris computer code ev8 [27l. I2H]. 

A. SCMF results 

The starting point of our study is presented in Fig. [2] where we show HF+BCS energy surfaces versus the axially 
symmetric mass quadrupole moment (Q) for the nuclei 20 Ne, 24 Mg and 36 Ar. These surfaces are calculated for the 
SLy4 force plus the zero-range density-dependent force (|27p by constraining the value of the mass quadrupole (Q). 
Results are shown for pairing strengths of g DFT = 900 and 1000 MeV-fm 3 . These energy surfaces are useful in finding 
the global minimum of the HF+BCS energy as well as providing the energy of the spherical solution. Unlike the 
HF surface of an open-shell nucleus, the (Q) — configuration is always a stationary point of the HF+BCS energy 
surface. 

In the HF approximation, the nucleus 20 Ne is prolate with a deformation energy of 5.6 MeV and (Q) m i — 84 
fm 2 [ID]. In HF+BCS, when including a pairing channel in the DFT force with a strength of <7 DFT = 900 MeV-fm 3 , 
it remains prolate but the deformation energy and quadrupole moment are reduced to E^+ BCS = 0.14 MeV and 
(Q)mf — 57.8 fm 2 , respectively. The deformed state is paired with ^^j^f = —3.3 MeV. However, i? pa i r is more 

negative in the spherical state E^^^ s — —6.16 MeV, lowering the mean-field energy of the spherical configuration 
more than in the deformed configuration, and hence the overall reduction of deformation energy. The second energy 
scale (|2~41 that is necessary for the construction of our map is given by <5-E^ a F r +BCS = 2.86 MeV. The results are 
tabulated in Table [Til 

As g DFT increases beyond 900 MeV-fm 3 , a shape transition occurs in 20 Ne to a spherical shape, and for g DFT = 1000 
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FIG. 2: HF+BCS energy curves for the nuclei 20 Ne (left column), 24 Mg (middle column) and 36 Ar (right column) as functions 
of the axially symmetric mass quadrupole moment (Q). The corresponding surfaces are shown for DFT pairing strengths of 
<? DFT = 900 and 1000 MeV-fm 3 . The SLy4 parametrization of the Skyrme interaction is used in the particle-hole channel while 
a zero-range density-dependent interaction is included in the pairing channel. All energies are measured with respect to the 
global minimum of the corresponding surface. 



TABLE II: HF+BCS results for the nuclei 20 Ne, 24 Mg and 36 Ar, using the SLy4 Skyrme interaction plus a zero-range a 
density-dependent force (|27jl in the pairing channel, are compared with HF+BCS results for the mapped Hamiltonian (|20|l . 
The coupling constants gQ and gp in (120 are determined by matching the deformation energy (|22|) and pairing energy excess 
flU) of the SCMF theory with those of the CISM. Energies are in MeV, (Q) is in fm 2 and g Q is in MeV/fm 4 . For the mapped 
Hamiltonian gp is in MeV while for the SCMF gp = <? DFT is in MeV-fm 3 . The last column shows the correlation energy (|25|l . 
calculated from the ground-state energy of the CISM Hamiltonian with the coupling parameters shown in the Table. 



Nucleus 


Interaction 


9Q 


gp 


r.HF + B('S 


rnHP+BCS 
pair.sph 


pHP+BCS 


^HF+BCS 
pair,mf 




r ^HF+BCS 


(Q) 


(Q)mf 




20 Ne 


SLy4 




900 


-157.91 


-6.16 


-158.05 


-3.30 


0.14 


2.86 


34.8 


57.8 






V 2 ' + pairing 


0.0954 


0.4516 


-35.58 


-5.00 


-35.72 


-2.14 


fit 


fit 


35.7 




3.4 


24 Mg 


SLy4 




900 


-192.61 


-5.84 


-195.93 





3.32 


5.84 


59.4 


111.9 






fo' 2 ' + pairing 


0.0343 


0.4930 


-85.12 


-5.84 


-88.44 





fit 


fit 


58.4 




4.7 


24 Mg 


SLy4 




1000 


-193.27 


-6.67 


-195.93 





2.66 


6.67 


59.4 


111.9 






+ pairing 


0.0332 


0.5729 


-86.02 


-6.67 


-88.68 





fit 


fit 


58.0 




5.1 


36 Ar 


SLy4 




900 


-305.02 


-2.14 


-305.43 





0.41 


2.14 


-43.4 


-74.7 






pairing 


0.0992 


0.3128 


-281.15 


-2.14 


-281.56 





fit 


fit 


-38.8 




2.0 


36 Ar 


SLy4 




1000 


-305.26 


-2.48 


-305.43 


-0.53 


0.17 


1.95 


-39.7 


-68.9 






+ pairing 


0.1017 


0.3275 


-283.15 


-1.95 


-283.32 





fit 


fit 


-35.4 




1.9 



MeV-fm 3 the ground-state configuration is spherical. In the present work, we discuss the mapping for nuclei whose 
mean-field solution is deformed, so the case g UFT = 1000 MeV-fm 3 is excluded for 20 Ne. 

The nucleus 24 Mg is also prolate in the HF approximation with a deformation energy of 8.73 MeV and (Q) m f — 111.9 
fm 2 [l8|. When a zero-range density-dependent force is included in the DFT pairing channel, 24 Mg remains prolate 
with the same deformation of (Q) m { = 111.9 fm 2 for both values of <? DFT . This is since the deformed HF+BCS ground 
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state is found to have n — 0. However, the spherical HF+BCS energy is lowered by the pairing force, decreasing the 
deformation energy to E^ { +BCS = 3.32 and 2.66 MeV for 5 DFT = 900 and 1000 MeV-fm 3 , respectively. In these 
cases, the pairing energy excess (|24|) is the absolute value of the spherical pairing energy and is observed to increase 
for larger values of the DFT pairing strength g DFT (see Table HI)) . 

Finally, the nucleus 36 Ar is oblate in the HF approximation with a deformation energy of 2.27 MeV and (<2)mf = 
—74.7 fm 2 (l8| . and it remains oblate with the same deformation in HF+BCS for ff DFT = 900 MeV-fm 3 (no pairing 
in the deformed solution). The pairing energy of the spherical configuration increases slightly with g DFT . However, 
the deformed solution becomes weakly paired for g BFT = 1000 MeV-fm 3 , and the ground-state quadrupole moment 
is reduced in magnitude to (Q) m f = —68.9 fm 2 . Consequently, we observe a small decrease in the value of SE^ ai ^ B 
(see Table HU. 

The DFT calculations provide us with the single-particle mean-field Hamiltonian of the deformed configuration, 
necessary for the construction of and in the effective CISM Hamiltonian ([20]) . as well as with the two energy 
scales E^ { +BCS and ^-Ep^^ 303 required for the determination of the two coupling parameters gQ and gp in (|20p . 

B. Truncation effects and operator rescaling 

The next step in constructing the map involves a transformation to the spherical basis and truncation of the 
spherical shell model space to the valence shells. This truncation was absent in our USD studies in Sec. IIII| since the 
original space was already truncated. Here the situation is different since the DFT includes a large number of spherical 
orbitals and it is necessary to truncate to the spherical sd shell. Thus, before we proceed with the construction of the 
CISM Hamiltonian, we study in this section truncation effects on the density matrix and pairing tensor, as well as on 
the quadrupole moment. 

The HF+BCS calculations determine the density matrix p and pairing tensor k (separately for protons and neu- 
trons) of the deformed mean-field ground state. We transform these matrices to the spherical DFT basis using Eqs. (|3]) 
and, as long as this spherical basis is complete, they continue to satisfy the consistency conditions (g]) and ([5]). 

Once the model space is truncated, the matrices p and k are projected on the valence sd shell to give PpP and 
PkP. We expect the mapping procedure to work well when the projected density matrix and pairing tensor satisfy 
approximately the consistency relations (fTB) and © . 

When the deformed HF+BCS solution is unpaired, k = and relations ([5|) reduce to the usual HF consistency 
relation p 2 = p. For such cases, this relation together with (|16p were already found to hold to a very good approxi- 
mation in Ref. 18J. Hence, we need to study truncation effects only in cases where the deformed mean-field solution 
is paired [k ^ 0). These cases are 20 Ne for 5 DFT = 900 MeV-fm 3 and 36 Ar for 5 DFT = 1000 MeV-fm 3 (see Table HJ. 

Rather than studying relations ((SJ), it is more convenient to study the equivalent requirement that the projected 
generalized density matrix TL is a projector in twice the number of dimensions [see Eq. I|14p]. Relation (fT4")) is satisfied 
exactly when 7Z has eigenvalues that are either 1 or and trace that is twice the number of valence single-particle 
orbitals. Table Hill lists the eigenvalues and trace of the matrix P7?.P. The dimension of PTZP in the sd shell is 24 for 
each type of nucleon. However, the eigenvalues of P7?.P are pairwise degenerate (because of time-reversal symmetry) 
and we list only the 12 distinct eigenvalues (for either protons or neutrons). We observe that all eigenvalues are very 
close to either 1 or 0, so that the projected generalized density matrix satisfies the HF+BCS consistency condition 
1Z 2 = TZ to a good accuracy. 



TABLE III: Eigenvalues and trace of the truncated generalized density matrix P7?.P of the deformed ground state in the SCMF 
theory of the SLy4 Skyrme energy functional plus a zero-range density-dependent force. The last column shows the trace of 
the truncated density matrix PpP. Results are shown for those cases in Table HT1 that have a paired deformed solution (k / 0). 



Nucleus 


Interaction 


,9 DFT 


nucleon 


Eigenvalues of P7?.P 




triPpP) 


20 Ne 


SLy4 


900 


neutron 
proton 


-6 x 10~ 6 -1 x 10" 6 4 x 10" 6 2.6 x 10" 5 1.7 x 10" 4 0.003 
0.986 0.991 0.991 0.992 0.995 0.996 

-3 x 10~ 6 1 x 10~ 6 5 x 10" 6 7.1 x 10" 5 3.8 x 10~ 4 0.003 
0.983 0.985 0.991 0.991 0.995 0.996 


23.824 
23.774 


1.986 
1.985 


36 Ar 


SLy4 


1000 


neutron 
proton 


-2 x 10~ 6 1.1 x 10" 5 4.1 x 10" 5 8.2 x 10~ 5 1.9 x 10" 4 8 x 10" 4 
0.993 0.995 0.996 0.997 0.998 0.999 

-1 x 10~ 5 -5 x 10~ 6 8 x 10~ 6 2.3 x 10 -5 8.5 x 10~ 4 
0.993 0.995 0.996 0.997 0.998 0.999 


23.918 
23.914 


9.964 
9.956 
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Another condition is (fT6| (i.e., the trace of the projected density PpP should be approximately equal to the valence 
number of nucleons of each type). When k = 0, this condition was already found to be well satisfied in Ref. [H. For 
cases with k ^ 0, the respective traces are tabulated in the last column of Table Mil We observe that condition 
(fTo| is satisfied to a good accuracy also when the ground-state solution is paired. We note that while the number of 
particles is fixed in the HF approximation, it is only the average number of particles that is fixed in the HF+BCS 
approximation. We also note that when the deformed ground state has /t/0, the eigenvalues of p correspond to the 
BCS parameters v\ and they can differ significantly from both 1 and 0. 

The SCMF theory allows us to determine the rescaling of one-body observables in the CISM model space by 
comparing their expectation values in the complete and truncated spaces. A good example is the mass quadrupole 
operator. In Table [TT1 we compare its expectation value (Q) m f in the complete space with its expectation value (Q) in 
the truncated space (sd shell). The ratio (Q)m{/(Q) determines the rescaling factor of the mass quadrupole moment 
when working in the truncated CISM space. In cases where k = 0, these scaling factors coincide with those found in 
the HF theory 18]. An interesting issue is the dependence of the rescaling factor on the strength <? DFT of the DFT 
zero-range density-dependent force, once the latter is sufficiently strong to induce pairing in the deformed ground- 
state solution. For 20 Ne, we find that the quadrupole operator should be rescaled by 1.66 for <7 DFT = 900 MeV-fm 3 , 
compared with 1.71 in the absence of a pairing channel in the DFT. This is a slight reduction of the rescaling factor 
by about 2.3%, even though the quadrupole moment itself is reduced by about 31%. Thus, the rescaling of the mass 
quadrupole operator in the truncated CISM space seems to be almost independent of g DFT . 



C. CISM Hamiltonian 

Once the monopole and quadrupole parts (h^ and hf^) of the projected mean-field Hamiltonian are determined, 
we can construct the effective CISM Hamiltonian (|2"0"1) except for the two unknown coupling constants gQ and gp. As 
already discussed in Sec. Ill CI these coupling constants are determined by matching simultaneously the DFT energy 
scales -Bdef +BCS and <^pah +BCS with similar energy scales for the CISM Hamiltonian when the CISM is solved in 
the HF+BCS approximation. The values of gQ and gp for the various cases studied here are listed in Table ITT1 In 
general, the spherical and deformed pairing energies of the mapped CISM Hamiltonian differ from the corresponding 
DFT pairing energies. For example, in 20 Ne with g DFT = 900 MeV-fm 3 , we find the values of E^+^ S = -5.00 

MeV and £" F r + ™f° S = -2.14 MeV in a mean-field solution of the CISM Hamiltonian, while the DFT values are 

JJcill ,1111 ' 

£pahtsph S = - 6 - 16 MeV and ^pa^rtrnf 053 = - 3 - 30 MeV ( see Table TO • u is onl y the difference ^p a F r +BCS in pairing 
energies between the deformed and spherical solutions that is matched (2.86 MeV in this case). 

We next study how well can the CISM Hamiltonian reproduce one-body observables such as the mass quadrupole. 
In Table [H] we compare the value of (Q) calculated when the CISM Hamiltonian (f2T))) is solved in the mean field (i.e., 
HF+BCS) approximation, with its SCMF value in the truncated space (sd shell). For 20 Ne at g DFT = 900 MeV-fm 3 , 
we find the CISM value of (Q) — 35.7 fm 2 as compared with the DFT value of (Q) — 34.8 fm 2 in the truncated sd 
space. We find a similarly close agreement for 24 Mg. The largest deviation (about 11%) is found for 36 Ar. 

Other relevant one-body quantities are the occupations (fij) of the spherical orbitals. In Table HVl we compare their 



TABLE IV: Neutron and proton spherical occupations (rij) for the nuclei 20 Ne, 24 Mg and 36 Ar in the sd-shell valence space. 
The DFT (SLy4 with a zero-range density-dependent force in the pairing channel) results are compared with the occupations 
of the mapped CISM Hamiltonian (|20[) both in the HF+BCS approximation and in shell model calculations. For the SLy4 
interaction gp = g DFT ■ 



Nucleus 


Interaction 


9Q 


gp 


n0d5/2 


nlsl/2 


n0d3/2 


p0d5/2 


plsl/2 


p0d3/2 


20 Ne 


SLy4 




900 


1.713 


0.221 


0.052 


1.709 


0.225 


0.052 




h {2) + pairing (HF+BCS) 


0.0954 


0.4516 


1.736 


0.216 


0.048 


1.731 


0.222 


0.048 




U 2) + pairing (CISM) 


0.0954 


0.4516 


1.730 


0.226 


0.045 


1.724 


0.230 


0.046 


24 Mg 


SLy4 




900 


3.293 


0.408 


0.215 


3.281 


0.401 


0.215 




U 2) + pairing (HF+BCS) 


0.0343 


0.4930 


3.489 


0.304 


0.207 


3.487 


0.306 


0.207 




h (2) + pairing (CISM) 


0.0343 


0.4930 


3.441 


0.299 


0.260 


3.435 


0.303 


0.263 


36 Ar 


SLy4 




1000 


5.934 


1.793 


2.237 


5.930 


1.786 


2.239 




h (2) + pairing (HF+BCS) 


0.1017 


0.3275 


5.995 


1.888 


2.117 


5.995 


1.884 


2.122 




h (2) +pairing (CISM) 


0.1017 


0.3275 


5.976 


1.836 


2.188 


5.975 


1.831 


2.194 
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FIG. 3: The occupations (hj) of the valence spherical orbitals as a function of the single-particle energy e - for the nuclei 
20 Ne, 24 Mg and 36 Ar (see Table HV)l . The DFT results (diamonds) are compared with the results in the HF+BCS ground-state 
solution of the CISM Hamiltonian (|20|l (triangles). Also shown are the occupations calculated in a shell model approach for 
the effective CISM Hamiltonian (squares). Open (solid) symbols are for protons (neutrons). 

DFT values (fij) = J2 m Pjm,jm with the respective spherical occupations of the CISM when solved in the HF+BCS 
approximation. We also list in Table |IV] the CISM occupations as calculated by the shell model code oxbash. Results 
are shown for 20 Ne and 24 Mg at 3 DFT = 900 MeV-fm 3 , and 36 Ar at 5 DFT = 1000 MeV-fm 3 . Even in cases where the 
deformed ground state is unpaired, the CISM occupations are expected to differ from their values in Ref. [H|. The 
reason is that the CISM Hamiltonian studied here contains an additional pairing term. The occupations computed in 
all three methods are in rather close agreement. These various spherical occupations are also shown in Fig. [3] versus 
the corresponding spherical single-particle energy . 

Finally, we compare the BCS occupations v\ of the deformed solution of the mapped CISM Hamiltonian (defined in 
the truncated space) with the v\ of the original SCMF ground state (defined in the complete space). Figure [4] makes 
this comparison for 20 Ne at g DFT = 900 MeV-fm 3 . The v\ of the mapped theory are slightly enhanced (suppressed) 




MMeV) 

FIG. 4: Neutron BCS occupations v\ of the deformed ground-state solution of the CISM Hamiltonian (|20[) (solid squares) 
are compared with the BCS occupations in the deformed basis of the DFT ground-state solution (open squares). These 
BCS occupations are plotted versus the deformed single-particle HF energies es. Results are shown for 20 Ne using the SLy4 
parametrization of the Skyrme force and a zero-range density-dependent force with <? DFT = 900 MeV-fm 3 . 



14 



below (above) the Fermi energy when compared with their DFT values. In the DFT, core orbitals below the sd shell 
have tijRil, while orbitals above the valence sd shell have v\ « 0. 

D. Correlation energies 

An important application of the DFT to CISM map is the calculation of correlation energies. We calculate the 
ground-state energy E gs of the mapped Hamiltonian (|20|) in the shell model approach and determine the correlation 
energy using Eq. (|25|) . The correlation energies for the three nuclei 20 Ne, 24 Mg and 36 Ar are tabulated in Table HT1 
In general, the correlation energy of the mapped Hamiltonian is found to increase with 5 DFT . In Fig. [5] we show the 
correlation energy versus mass number A for <? DFT = 900 MeV fm 3 . We compare our results (solid squares) with the 
correlation energies obtained in Ref. Jjl where the effective interaction includes only a quadrupolar component (open 
squares). The inclusion of pairing correlations provides an additional gain in correlation energy and the observed 
trend versus mass number A is similar to the trend observed in the USD results shown in Fig. Q] 




20 24 28 32 36 
A 

FIG. 5: Correlation energies -E C orr [see Eq. 125)) ] versus mass number A for the N = Z sd-shell nuclei 20 Ne, 24 Mg and 36 Ar, 
calculated using the mapped CISM Hamiltonian (|20p (solid squares). The SCMF theory used to construct the map is the SLy4 
Skyrme force plus a zero-range density-dependent force in the pairing channel of strength <? DFT = 900 MeV fm 3 . The open 
squares are the correlation energies found from a mapped CISM Hamiltonian that includes quadrupolar correlations alone 

Our correlation energies of 5.1 MeV and 1.9 MeV for the nuclei 24 Mg and 36 Ar, calculated for g DFT = 1000 MeV 
fm 3 , are in very close agreement with the values of 5.1 MeV and 2.0 MeV obtained in Ref. [TH starting from the same 
DFT but using an angular-momentum projected GCM. The correlation energy of 3.4 MeV we find for 20 Ne (using a 
somewhat reduced value of <? DFT ) is also consistent with the value of 3.1 MeV reported in Ref. [3. 

V. CONCLUSION 

In this work we have mapped an SCMF theory onto an effective CISM Hamiltonian in the presence of both 
quadrupolar and pairing correlations. When compared with Ref. I la. this is an important step forward in constructing 
a general map that would take advantage of the global validity of the energy density functional and the higher accuracy 
of the CISM approach. The main differences compared with the study initiated in Ref. [l8[ are: (i) The SCMF is now 
based on an HF+BCS approximation (rather than HF), and the ground-state consistency conditions are recast in 
terms of both the density matrix and the pairing tensor, (ii) The effective Hamiltonian depends now on two coupling 
strengths, and consequently two energy scales must be matched to determine these coupling parameters. One energy 
scale is the deformation energy, similar to the one used previously but now redefined in terms of the HF+BCS energies. 
The second scale is determined by the difference of average pairing energy in the deformed ground state and spherical 
configurations. The new map constructed here provides improved correlation energies that are consistent with other 
methods that include quadrupolar and pairing correlations beyond the mean field. 
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The DFT application discussed here is based on a Skyrme force in the particle-hole channel plus a zero-range 
density-dependent interaction in the pairing channel. One could alternatively start from a mean-field calculation with 
a finite-range Gogny force. Such an interaction includes pairing correlations consistently and is usually treated in the 
Hartree-Fock-Bogolyubov (HFB) method. The map constructed here can be easily generalized to the Gogny DFT by 
replacing the HF+BCS approximation by the HFB approximation, both when solving the DFT and in the mean-field 
solution of the CISM Hamiltonian. 



As already discussed in Ref. [18[ much remains to be done. Nuclei that have no spherical pairing solution can 

be treated by introducing a constraining pairing field. The addition of an interaction term built from an isovector 
quadrupole operator would enable us to generalize our studies to N ^ Z nuclei and carry out systematic studies. A 
major challenge is the application to the heavy deformed nuclei. Although the method presented here can in principle 
be adapted to such nuclei, outstanding issues are the ability of a spherical shell model to describe nuclei with large 
deformation in a truncated space and the practical issue of solving the CISM Hamiltonian in large model spaces. 
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